A recursive method for the density of states in one dimension 
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We derive a powerful yet simple method for analyzing the local density of states in gapless 
one dimensional fermionic systems, including extensions such as momentum dependent interaction 
parameters and hard-wall boundaries. We study the crossover of the local DOS from individual 
density waves to the well-known asymptotic powerlaws and identify characteristic signs of spin charge 
separation in possible STM experiments. For semi-infinite systems a closed analytic expression is 
found in terms of hypergeometric functions. 
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The density of states (DOS) is a central quantity in 
the study of electronic condensed matter systems. The 
corresponding expression for the local DOS 

p(u),x) = ^2\{u m \&(x)\0)\ 2 5(u-u m ) 
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describes the probability of inserting an electron at a 
given energy oj and can be probed in tunneling exper- 
iments or, averaged over a range x, in photoemission ex- 
periments. In two or more dimensions the DOS is typ- 
ically peaked at renormalized single particle excitation 
energies. In one dimensional systems, however, interac- 
tion effects are enhanced, so that the DOS is determined 
by collective many-body states instead Accordingly, 
there has been considerable theoretical interest in analyz- 
ing the one dimensional DOS for several decades [lrtZ|. 
Typical features that have been predicted are separate 
spin and charge density excitations and a characteristic 
depletion of the DOS at low energies and near bound- 
aries, which have been seen experimentally in some spe- 
cial cases, such as cleaved edge overgrowth wires su- 
perstructures on surfaces || , and carbon nanotubes flCj| . 
Recent theoretical activities in the field have produced 
notable advances in the areas of non-linear corrections 
, numerical simulations 12|, |13[ , and applications of 
exact methods 

On the other hand, there remain a number of open 
questions, especially in regards to the applicability of 
typical effective low energy theories such as the Luttingcr 
liquid formalism to more realistic systems and models. In 
fact a number of energy and length scales can affect the 
behavior, including perturbations from other degrees of 
freedom, longer range electron-electron interactions, im- 
purities, and finite system sizes, so that realistic systems 
are never truly scale invariant. The renormalization due 
to higher order perturbing operators leads to an energy 
and momentum dependence of the interaction parame- 
ters, so that a description in terms of a single constant 
Luttingcr liquid parameter is in general not adequate. 



Accordingly, the low energy theory is altered substan- 
tially and the generically predicted powerlaw behavior 
with energy, momentum, and position is changed into a 
more complicated behavior of the DOS. 

In this paper, we address the question of how to gen- 
erally calculate the local DOS including complications 
which make the central interaction parameters effectively 
energy and momentum dependent. The calculation is 
based on the Fourier transform in Eq. (JlJ for vertex cor- 
relation functions in finite systems, which yields an ex- 
pansion in delta-functions for the discrete DOS with coef- 
ficients that follow a recursion relation for arbitrary mo- 
mentum dependent interaction parameters. It is shown 
that boundaries cause a natural scale dependence in the 
description of the local DOS, which leads to the crossover 
from boundary to bulk behavior. Using a continuous de- 
scription, we obtain a closed analytic expression of the 
local DOS as a function of energy and position. In or- 
der to identify spin- and charge separation by scanning 
tunneling microscopy (STM) it is useful to analyze the 
spatial Fourier transform of the local DOS. 

The starting point is the general expression of chiral 
fermionic fields in terms of vertex operators in the normal 
ordered form 



0*(x,t) :-- 



where Ai are linear combinations of bosonic annihilation 
operators, e.g. of the form Ag(x) — ate lk ' x bf+fiie~ lkix b^ 
for a periodic system with length L, energy spacing 
Aw = and k t = ^ 0. Such operators are 
used to represent left- or right moving fermion operators 
ip\/ R - In the case of several fermion channels (e.g. spin 
and charge) the chiral fermion field will be the product 
of operators 0^(x,t) for each channel separately, all of 
which have the form in Eq. ^ , but with different energy 
spacings Alo. The magnitude of the prefactor c is typi- 
cally unknown unless a comparison with exact results can 
be made. Zero mode terms have been omitted in Eq. ([2]) 
since they only shift the spectrum of the DOS. 

The form of the vertex operators in Eq. ([2]) is believed 
to apply to gapless interacting fermion systems in one 
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dimension. However, it is important to note that the 
expression (|2|) already implies that the possible energy 
levels are assumed to be evenly spaced at cj = £Aoj rela- 
tive to the Fermi energy ep, which is generally not exactly 
justified in realistic systems and is also explicitly violated 
in any finite lattice model. The reason why bosonization 
remains useful is that the central information about the 
electron-electron interactions is encoded in the Bogoli- 
ubov rotations, i.e. the exact form of the linear combina- 
tions of bosons Ai(x) in Eq. ([2]). These do not crucially 
depend on the assumption of an equally spaced spectrum, 
since they are derived in momentum space. Therefore, 
the typical bosonization theory does not describe the ex- 
act energy locations of the levels very accurately, but it 
appears to predict the spectral weights surprisingly well, 
as can also be seen by numerical studies |13| ■ 

The general time correlation function can be calculated 
from Eq. © 
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(O(x,t)O*(x,0)) = \c\ 2 expC£je- iiAut le(x)), (3) 
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where ji(x) = [A,(x),A'f(x)]. In accordance with the pe- 
riodicity in time, the Fourier transform in Eq. (|T|) yields 
an expansion in delta-functions for the DOS 
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By partial integration a recursion formula for the indi- 
vidual spectral weights p m is obtained 
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with po = |c| 2 . This simple but central result is the main 
instrument for analyzing the DOS in what follows. Be- 
sides the obvious simplicity of calculating the DOS with 
Eq. (0, it has the major advantage of being free from 
divergences and regularization schemes. Moreover, it is 
possible to analyze a dependence of the j£ on the prin- 
ciple energy quantum number I which will be especially 
useful for boundaries and long range interactions. 

Let us first re-examine the simplest example of a Lut- 
tinger liquid, which is a short-range interacting single 
channel system with periodic boundary conditions. The 
commutator 7 = [At, A' e ] = h + K) is then given in 
terms of the Luttinger parameter K independent of t. In 
this case, Eq. ([5]) is solved analytically by 
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This is a well known result [3[, which will be useful 
later. However, many states can potentially contribute to 
each spectral weight p m , the number of which increases 
with the partitions of to. As mentioned above, those 



states are in general not exactly degenerate at the energy 
toAw. Therefore, the weights p m are spread in energy 
[ill, which is also the case in all examples that follow. 
However, the total value of p m is expected to be rather 
accurate 

Next, we consider fermions in a finite system with 
hard-wall boundary conditions ^(0) = ^(L) = 0. In this 
case, the operators in Eq. © are of the form Ai(x) = 



(a e e lk ' x + Pie- lkeX )bf with k e = ^ [17J. The Green's 
function splits into a uniform and a 2kpx oscillating part 
(*(x,t)&(x,0)) = (G R (x,x,t) - e i2k ^G R (x,-x,t)) + 
c.c, where G R (x,y, t) = {^ R (x,t)^ R (y,0)) B®M- For 
spinful fermions those correlation functions have the form 
of Eq. © for the spin factor (K s , v s ) and the charge fac- 
tor (K c , v c ) separately 



G R {x,x,t) = \c x \ 2 ] [ exp( 



E 



1 



7"M) (7) 



G R (x,-x,t) = \c x \ 2 ] [ cxp(W i; 
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7°T(*)), 



where the commutators are dependent on position x 
and mode I in this case 7 "" 1 (a;) = + ^f- + 

(ik - ^)cos(2^x) and ~/™(x) = - ^) ' 
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-^J cos(2fc^x) -I- I sin(2fc^.T) with Aw„ = 

for v = s,c. Also the prefactor is position dependent 

\c x \ 2 oc Iliy( sm 17) from normal ordering. Us- 
ing Eq. (O for the uniform/oscillating and spin/charge 

, n uni/osc 1 \-^m uni/osc uni/osc ., . 

parts separately p„, m ' = — }_^ e=1 p v ^_ t J„ tt , it is 
possible to solve the Fourier transform ([1]) even for this 
more complicated case. The local DOS at the energy 
lu = m c Auj c + m s Acj g is then simply the folded product 
of spin and charge, i.e. 

Pu/o= E P^ osc P^Hu-^cAiOc-msAus) (8) 

m c .m s 

and p(uj, x) = (p u (w, x) — e l2kpX p (ui, x)) + c.c. 

Equations §5§ and |5]) can be evaluated very quickly 
with a few lines of code for a large number of levels and 
different parameters. The analytic expressions for each 
level are superpositions of density waves up to wavenum- 
bers corresponding to m s and m c . The level density in- 
creases with uj, according to the number of ways of choos- 
ing m s and to c for a given energy interval in Eq. ([8]) . It is 
therefore useful to average over a small energy window, 
which results in a typical local DOS as a function of x 
and w as shown in Fig. [1] for K c = 0.531, K s = 0.618 
and v c /v s = 1.618, where the oscillating part has been 
omitted. The solid lines in Fig. [T] show the characteris- 
tic standing density waves with higher wavenumbers as 
the energy increases, with scaling of features along hy- 
perbolas with uj = mv c J x, while v s plays a lesser role. In 



FIG. 1: (Color online) p u (uj,x) in units of \c x= l/2\ 2 for K c — 
0.531, K s = 0.618 and v c /v s — 1.618 averaged over an energy 
interval of 0.3 in units of Acj c . 



FIG. 2: (Color online) Spatial Fourier transform of the local 
DOS with parameters as in Fig. [T]in arbitrary units. Inset: 
Cut at fixed energy. 



contrast, for non-interacting systems K c = K s = 1 and 
Fig. Q] is completely flat. 

In order to get clear evidence of the spin and charge 
separation it is necessary to consider the oscillating part, 
which is best analyzed by a Fourier transform as first 
considered in Ref. 0. In Fig. H we have plotted the 
absolute value of the Fourier transform in space of the 
local DOS including uniform and oscillating parts. It is 
possible to identify separately dispersing spin and charge 
features starting at 2k p. In addition there is a third 
maximum at 2&f for all energies. The uniform DOS at 
k = shows the typical powerlaw increase. For spin- 
independent interactions K s = 1,K C < 1 the results for 
Figs. [T] and [5] look qualitatively very similar (not shown), 
except for a missing slowly dispersing maxima which can 
barely be detected in Fig. [2] starting from k = 0. 

A further analytical analysis of the local DOS is possi- 
ble in the continuous limit of the level spacing xm v 3> L. 
In this case, Eq. ([5]) becomes an integral equation 

/>LU 

ojp(ui,x)= / p((*)', x) 7(cj — a/, x) dui' , (9) 
Jo 

which again holds for the oscillating/uniform and 
spin/charge parts separately, where 7(w, x) are the corre- 
sponding continuous functions. Let us focus on the uni- 
form charge part close to a single edge, i.e. 7(w, x) = 

ik + if + (ife - t 1 ) cos (^) from above - After 

rewriting Eq. ([9]) in terms of the dimensionless variable 
y = ujx/v c , differentiating three times, and using the fact 
that 7" = 4( j^- + — 7), it is possible to obtain a third 
order differential equation, which in turn can be solved 



by the following hypcrgcomctric function 
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Y^oo r(a+fc) 2* 

r(a) ' 2^,k=o r(6i+fc)r(h 2 +fe) k\ • 
The dependence on position is not yet specified. The full 
local DOS from several channels can then be obtained 
by folding according to the continuous version of Eq. (|8]) 
which for K s = 1 again results in a hypergeometric func- 
tion 



Pu = 



f{x) v t f pT\v')pT - y')) dy' (11) 
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as shown in Fig. [3l The x-dependence f(x) is fixed by 
the asymptotic behavior near x = which yields the 
expression ((6|) with 7 = + i . Therefore an exact an- 
alytic formula (fT2"j) has been derived without the need to 
evaluate the Fourier transform of correlation functions in 
Eq. ([T]), which involves a complicated contour integration 
and has so far not been possible analytically. Figure 
[3] also shows the DOS for all m s and m c quantum num- 
bers from the discrete recursion formula ([S]). Clearly, the 
individual weights drop with increasing w, but the aver- 
aged behavior follows the analytical prediction (fTS]) . The 
boundary and bulk exponents of the effective powcrlaws 



p cx Lo a are given by a ~ ^ v 7^ 
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1 expanded for 



small x or averaged over large x, respectively, which are 
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FIG. 3: (Color online) The uniform spectral weights p u in 
units of \c x \ 2 near a boundary as a function of ljx/v c for 
K c = 0.531, K a = 1, and v c /v a — 1.618. The large number of 
spin and charge levels for x — 0.1L (blue dots) follow Eq. (|T2j) 
(red line) after averaging. 



however only rough approximations compared to the an- 
alytic expressions above. 

So far, we have considered a constant Luttinger pa- 
rameter K, where the scale dependence of the DOS stems 
only from the boundary. However, for longer range inter- 
actions, Ki itself is expected to be dependent on £ even in 
the bulk, e.g Kg = 1/^/Vb hi | for the case of a Coulomb 
potential [18|. There has been a controversy in the lit- 
erature about the behavior of the corresponding DOS, 
for which an effective exponent p oc ui a ^ u} ' was postu- 



19. 20] 



latcd in Ref . [19( , which in turn lead to a debate 

Our analysis with = \ + KA using Eq. ([5]) now 

shows that over a limited range the true behavior can 
fit to the form postulated in Ref. 



19j by some choice 



of parameters, but that there is an additional unknown 
scale dependence as the authors also mention 19]. There- 
fore, the description in terms of a(to) is not complete, 
but phenomenologically useful. An analysis with help 
of Eq. (|5|) shows that corrections come from derivatives 
in which however cannot be summed exactly. Using 
the recursion formula it is possible to consider interac- 
tions Ki of arbitrary form with relative ease. Also for 
harmonically trapped fermions a description of a mode 
dependent parameter Kg has been predicted, albeit with 
a yet unknown dependence on £ [2l| . 

This brings us to an important application of the re- 
cursion formula for systems which are not perfect Lut- 
tinger liquids. As outlined in the introduction, poten- 
tially all realistic models contain some scale-dependence 
of the parameter Kg, which is however not a priori known. 
Therefore, it will be very beneficial to use the recursion 



formulas ([5]) in reverse: Knowing the lowest m spectral 
weights pi, it is easily possible to uniquely determine 
the first m — 1 interaction parameters "/£. The analo- 
gous statement is also true for continuous spectra. This 
is especially promising considering the fast progress in 
numerical [II, [Tj|, exact [14j . and experimental (8HTfl] 
techniques in determining the DOS. 

In summary, we have obtained a straight-forward tool 
for analyzing the local DOS in one dimensional inter- 
acting systems. The crossover from individual levels to a 
continuous spectrum has been studied in detail for future 
comparison with possible STM experiments. A closed an- 
alytic expression was obtained for electrons in a system 
with hard-wall boundaries. Long range interactions can 
also be considered. The inversion of the recursion for- 
mula is useful to calculate the interaction parameters di- 
rectly from the low energy spectral weights, which in turn 
can be obtained from numerical, exact, or experimental 
methods. In this way the theory can be generalized to 
systems which are not perfect Luttinger liquids, such as 
harmonically trapped ultra-cold fermions. 
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